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Abstract 


The centrality dependence of the charged-particle pseudorapidity density measured with ALICE in 
Pb-Pb collisions at = 2.76 TeV over a broad pseudorapidity range is presented. This Letter 

extends the previous results reported by ALICE to more peripheral collisions. No strong change of 
the overall shape of charged-particle pseudorapidity density distributions with centrality is observed, 
and when normalised to the number of participating nucleons in the collisions, the evolution over 
pseudorapidity with centrality is likewise small. The broad pseudorapidity range (—3.5 < rj < 5) 
allows precise estimates of the total number of produced charged particles which we find to range 
from 162±22(syst.) to 17170±770(syst.) in 80-90% and 0-5% central collisions, respectively. The 
total charged-particle multiplicity is seen to approximately scale with the number of participating 
nucleons in the collision. This suggests that hard conttibutions to the charged-particle multiplicity 
are limited. The results are compared to models which describe dVch/drj at mid-rapidity in the 
most central Pb-Pb collisions and it is found that these models do not capture all features of the 
distributions. 
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1 Introduction 

The measurement of the charged-particle pseudorapidity (rj) density distribution in heavy-ion collisions 
provides insight into the dominant particle production mechanisms, such as parton fragmentation [1] and 
the observed phenomenon of limiting fragmentation [2]. The unique capability of ALICE to perform 
such measurements from large to small overlaps of the colliding nuclei over a broad pseudorapidity range 
allows for significant additional information to be extracted e.g., the total number of charged particles 
and the evolution of the distributions with centrality. 

The charged-particle pseudorapidity density (dAch/dq) per se does not provide immediate understanding 
of the particle production mechanism, but as a benchmark tool for comparing models it is indispensable. 
Various models [3-5] make different assumptions on how particles are produced in heavy-ion collisions 
resulting in very different charged-particle pseudorapidity density distributions — both in terms of scale 
and shape. Models may, for example, incorporate different schemes for the hadronisation of the produced 
quarks and gluons which leads to very different pseudorapidity distributions of the charged particles. 

The ALICE collaboration has previously reported results on the charged-particle pseudorapidity density 
in the 0-30% most central Pb-Pb collisions at = 2.76TeV over a wide pseudorapidity range [6], 
and in the 80% most central collisions at mid-rapidity (rj ss 0) only [7]. The ATLAS collaboration has 
reported on the charged-particle pseudorapidity density in the 80% most central events in a limited pseu¬ 
dorapidity range of |t] | < 2 [8]. Similarly, the CMS collaboration has reported on the same measurements 
in the 90% most central events at T] « 0, and for selected centralities up to ItjI < 2 [9]. 

In this Letter we present the primary charged-particle pseudorapidity density dependence on the event 
centrality from mid-central (30-40%) to peripheral (80-90%) collisions over a broad pseudorapidity 
range to complement results previously reported by ALICE in the 0-30% centrality range. Unlike pre¬ 
vious [6] , in the forward regions where the signal is dominated by secondary particles produced in the 
surrounding material, we use a data-driven correction to extract the primary charged-particle density. 

Primary charged particles are defined as prompf charged particles produced in fhe collision, including 
fheir decay producfs, buf excluding producfs of weak decays of muons and light flavour hadrons. Second¬ 
ary charged parficles are all ofher particles observed in fhe experimenf e.g., parficles produced fhrough 
inferacfions wifh maferial and producfs of weak decays. 

In fhe following secfion, fhe experimenfal sef-up will be briefly described. Section 3 ouflines analysis 
procedures and describes a dafa-driven mefhod fo isolate fhe number of primary charged parficles from 
fhe secondary parficle background af large pseudorapidify. Sysfemafic uncerfainfies are discussed in 
Sect. 4. In Sect. 5, the resultant charged-particle pseudorapidity density distributions are presented along 
with their evolution with centrality. Lurthermore we extract from the measured dVch/dt] distributions 
the total number of charged particles as a function of the number of participating nucleons. We finally 
compare the measured charged-particle pseudorapidity density to a number of model predictions before 
concluding in Sect. 6. 

2 Experimental setup 

A detailed description of ALICE can be found elsewhere [10, 11]. In the following we briefly describe 
the detectors relevant to this analysis. 

The Silicon Pixel Detector (SPD) is the inner-most detector of ALICE. The SPD consists of two cyl¬ 
indrical layers of 9.8 x 10® silicon-pixels possessing binary read-out. It provides a measurement of 
charged particles over | T] | < 2 using so-called tracklets — a combination of hits on each of the two lay¬ 
ers (1 and 2) consistent with a track originating from the interaction point. Possible combinations of 
hits not consistent with primary particles can be removed from the analysis, with only a small (a few %) 
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residual correction for secondary particles derived from simulations. The SPD also provides a measure¬ 
ment, by combining hits on its two layers, of the offset with respect to the interaction point, where the 
collisions occurred. IP = (0,0,0) is at the centre of the ALICE coordinate system, and IP^ is the offset 
along the beam axis. Einally, a hardware logical or of hits in each of the two layers provides a trigger for 
ALICE. 

The Eorward Multiplicity Detector (EMD) is a silicon strip detector with 51200 individual read-out chan¬ 
nels recording the energy deposited by particles traversing the detector. It consists of three sub-detectors 
EMDl, 2, and 3, placed approximately 320cm, 79 cm and —69cm along the beam line, respectively. 
EMDl consists of one inner type ring (li), while both EMD2 and 3 consist of inner (2i,3i) and outer type 
rings (2o,3o). The rings have almost full coverage in azimuth {(p), and high granularity in the radial ( 77 ) 
direction (see Tab. 1). 

The VO is the most forward of the three detectors used in this analysis. It consists of two sub-detectors: 
VO-A and VO-C placed at approximately 333cm and —90cm along the beam line, respectively. Each of 
the sub-detectors are made up of scintillator tiles with a high timing resolution. While the VO provides 
pulse-height measurements, the energy-loss resolution is not fine enough to do an independent charged 
particle measurement. In previous measurements, using so-called satellite-main collisions (see Sect. 3), 
one could match the VO amplitude to the SPD measurements to obtain a relative measurement of the 
number of charged particles. However, for collisions at jlP^I < 15 cm no such matching is possible, and 
the VO is therefore not used to provide a measurement of the number of charged particles in this analysis. 
The detector is used, in an inclusive logical or with the SPD, for triggering ALICE and to provide a 
measure of the event centrality [7]. 

Details on the coverage, resolution, and segmentation of the three used detectors are given in Tab. 1 . 


Detector 

dr(p 

5z 

Tj range 

SPDl 

12 pm 

100 pm 

—2.0 to 2.0 

2 

12 pm 

100 pm 

-1.4 to 1.4 

Detector 

A(p 

Ar 

J] range 

FMDli 

18° 

254 pm 

3.7 to 5.0 

2i 

18° 

254 pm 

2.3 to 3.7 

2o 

9° 

508 pm 

1.7 to 2.3 

3o 

9° 

508 pm 

-2.3 to-1.7 

3i 

18° 

254 pm 

-3.4 to -2.0 

VO-A 

45° 

34 to 186 mm 

2.8 to 5.1 

-C 

45° 

26 to 127 mm 

-3.7 to-1.7 


Table 1: Overview of the resolution (5), segmentation (A), and coverage of the detectors used in the analysis. The 
‘A’ side corresponds to z > 0, while the ‘C’ side corresponds to z < 0. The 7] range is specified for collisions with 
IP, = 0. 

3 Data sample and analysis method 

The results presented in this paper are based on Pb-Pb collision data at = 2.76TeV taken by 
ALICE in 2010. About 100000 events with a minimum bias trigger requirement [7] were analysed 
in the centrality range from 0% to 90%. The data was collected over roughly 30 minutes where the 
experimental conditions did not change. 

The standard ALICE event selection [12] and centrality estimator based on the VO-amplitude are used 
in this analysis [13]. We include here the 80-90% centrality class which was not present in the previous 
results [7]. As discussed elsewhere [13], the 90-100% centrality class has substantial contributions from 
QED processes and is therefore not included in this Letter. 

Results in the mid-rapidity region (| T] | < 2) are obtained from a tracklet analysis using the two layers of 
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the SPD as mentioned in Sect. 2. The analysis method and data used are identical to what has previously 
been presented [ 6 , 7]. 

The measurements in the forward region (|t]| > 2) are provided by the FMD. The FMD records the full 
energy deposition of charged particles that impinge on the detector. Since all charged particles that hit 
the FMD are boosted in the laboratory frame, the detection efficiency is close to 100% for all momenta. 
As reported earlier [ 6 ], the main challenge in measuring the number of charged primary particles in this 
region, is the large background of secondary particles produced in the surrounding material. Due to the 
complexity and the limited knowledge of the material distribution of support structures away from the 
central barrel, it has not been possible to adequately describe (on the few %-level) the generation of 
secondary particles in the forward directions within the precision of the current simulation of the AFICE 
apparatus. 

A suitable means to extract the number of primary charged particles was found by utilising collisions 
between so-called ‘satellite’ bunches and main bunches offset in intervals of 37.5cm along the beam¬ 
line. Satellite bunches are caused by the so-called debunching effect [14]. A small fraction of the 
beam can be captured in unwanted RF buckets, due the way beams are injected into the accelerator, and 
create these satellite bunches spaced by 2.5 ns. Collisions between satellite and main bunches can cause 
instabilities in the beam, and the FHC has taken steps to reduce the number of these kinds of collisions. 
AFICE has therefore not recorded collisions between satellite and main bunches before or after the Pb- 
Pb run of 2010. In satellite-main collisions the background of secondary particles was much smaller and 
much better understood since significantly less detector material shadows the forward detectors [ 6 ] . 

A study utilising these satellite-main collisions led to the publication of the measurement of the charged- 
particle pseudorapidity density in the 30% most central events over |t]| <5 [6 ]. The study was limited 
in centrality reach by the need to use the Zero-Degree Calorimeter (ZDC) for the centrality estimation 
for collisions between satellite and main bunches. The ZDC measures the energy of spectator (non¬ 
interacting) nucleons with two components: one measures protons and the other measures neutrons. The 
ZDC was located at about 114 m from the interaction point on either side of the experiment [10], and 
was therefore ideally suited for that study. The centrality determination capability of the ZDC is however 
limited to the 30% most central collisions [13]. 

For centralities larger than 30% the VO amplitude is used as the centrality estimator, which is available 
only for collisions at jlP^I < 15 cm — the so-called nominal interaction point corresponding to main 
bunches of one beam colliding with main bunches of the other beam. 

To extend the centrality reach of the dVch/dt] measurement, a data-driven correction for the number of 
secondaries impinging on the FMD has been implemented. For each centrality class C, we form the ratio 

d^ch/dt] I c.primary.satellite 

That is, the ratio of the measured inclusive charged-particle density from main-main collisions (lIP^I < 
10 cm) provided by the FMD to the primary charged-particle density from satellite-main collisions [ 6 ]. 
Here, ‘inclusive’ denotes primary and secondary charged particles i.e., no correction was applied to 
account for secondary particles impinging on the FMD. 

Note, that the correction is formed bin-by-bin in pseudorapidity, so that the pseudorapidity is the same 
for both the numerator and denominator. However, the numerator and denominator differ in the offset 
along the beam line of origin of the measured particles: For the numerator the origin lies within the 
nominal interaction region, while for the denominator the origin was offset by multiples of 37.5 cm. 

This ratio is obtained separately for all previously published centrality classes: 0-5%, 5-10%, 10-20% 
and 20-30%. The variation of Ec for different centralities is small (< 1%, much smaller than the precision 
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of the measurements). The weighted average 


E{^) 


IcACEc(t7) 

IcAC 


is used as a global correction to obtain the primary charged-particle pseudorapidity density 


d^ch 


dr] 


X,primary 


1 (Wch 

E(r]) dr] 


X ,inclusive,nominal 


where X stands for an event selection e.g., a centrality range. 


( 2 ) 

(3) 


The simulation-based correction S{r]) for secondary particles to the charged-particle pseudorapidity 
density in the forward directions is given by 


siv) 


Vinclusive,FMD(t] ) 
Vprimary.generated (t] ) 


(4) 


where Vinciusive,FMD is the number of primary and secondary particles impinging on the FMD — as 
given by the track propagation of the simulation, and Apnmaiy.generated is the number of generated primary 
particles at a given pseudorapidity. Complete detector-simulation studies show that three effects can 
contribute to the generation of secondaries, and hence the value of S(t]). These three effects are: material 
in which secondaries are produced, the transverse momentum (pj) distribution and particle composition 
of the generated particles, and lastly the total number of produced particles. Of these three the material 
is by far the dominant effect, while the pj and particle composition only effects S{r]) on the few percent 
level. The total number of generated particles has a negligible effect on S{r]). That is, the material 
surrounding the detectors amplifies the primary-particle signal through particle production by a constant 
factor that first and foremost depends on the amount of material itself, and only secondarily on the pj 
and particle composition of the generated primary particles. 


To estimate how much Ec{r]) itself would have changed if another system or centrality range was used to 
calculate the correction, S{r]) is analysed from simulations with various collision systems and energies. 
We find fhaf, even for large variations in parficle composifion and pj disfribufions, S{r]) only varies 
by up fo 5%. Re-weighling fhe particle composition and pj disfribufions from fhe various sysfems 
fo mafch produces consisfenf values of S{r]) ensuring fhaf fhe 5% variafions found were only due fo 
particle composifion and pj disfribufions differences. This uncerfainfy is applied to E(r]) fo accounf for 
all reasonable variations of fhe parficle composifion and pj disfribufions, which cannof be measured in 
fhe forward regions of ALICE. 


Figure 1 shows fhe comparison of fhe dafa driven correction E(r]) fo fhe simulafion-based correction 
S(t]) from Pythia [15] (pp) and a parameferisafion of fhe available ALICE resulfs [16, 17] for Pb- 
Pb collisions. The simulafed collisions are for fwo disfincf sysfems and span over almosf an order of 
magnifude in collision energy. The fofal number of produced particles in fhese simulafions span five 
orders of magnifude, and no dependence of S(r]) on charged-particle mulfiplicify is observed. 


By comparing ^(t]) fo S(t]) from simulations, one finds a good correspondence befween fhe fwo correc¬ 
tions except in regions where fhe maferial description in fhe simulafions is known fo be inadequafe. This, 
fogefher wifh fhe facf fhaf fhe numerafor and denominafor of Eq. 1 measure fhe same physical process, 
buf differ foremosl in fhe maferial fraversed by fhe primary particles, and hence fhe number of secondary 
particles observed, implies fhaf fhe correction E{r]) is universal. Thai is, Eq. 3 is applicable for any 
even! selection X in any collision sysfem or af any collision energy, where fhe produced mulfiplicify, px 
disfribufions, and particle composition is close fo fhe range of fhe simulated sysfems used fo sludy S{'r]). 


Nofe, for fhe previously published resulfs [6], which used salellife-main collisions, fhe simulafion-based 
approach for correcling for secondary particles i.e., applying S{'r]) direcfly, was valid. As mentioned 
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Figure 1: [Colour online] Comparison of data-driven to simulation-based corrections for secondary particles 
impinging on the FMD. Different markers correspond to different collision systems and energies, and the colours 
indicate the hve FMD rings. S{rj) is shown for Ocm < IPj < 2cm as an example, while E{rj) is independent 
of IPj (see also text). PYTHIA was used for pp collisions, and the Pb-Pb points are from simulation with a 
parameterisation which include the available ALICE data on particle composition and pj distributions. Black 
circles correspond to E{rj). 

above, in satellite-main collisions, the particles that impinge on the FMD traverse far less and better 
described material in the simulation of the AEICE apparatus. The use of a simulation-based correc¬ 
tion for secondary particles was in that analysis cross-checked by comparing to and combining with 
measurements from the VO and SPD [6]. Despite concerted efforts to improve the simulations by the 
Collaboration it has not been possible to achieve the same accuracy in ^(i]) for main-main collisions. 

Finally, the effect of variation of the location of the primary interaction point on E{ri) was studied. It was 
found, that the effect is negligible, given that the distribution of IP^ are similar between the numerator of 
Eq. 1 and right-hand side of Eq. 3, as was the case in this analysis. 

The method used in this analysis to extract the inclusive number of charged particles from the FMD is 
the same as for previous published results [6], except that the data-driven correction E{r\) — rather than 
a simulation-based one 5(1]) — is used to correct for secondary particles. 

4 Systematic uncertainties 

Table 2 summarises the systematic uncertainties of this analysis. The common systematic uncertainty 
from the centrality selection is correlated across T] and detailed elsewhere [13]. 

For the SPD measurements, the systematic uncertainties are the same as for the previously published 
mid-rapidity result [7], except for a contribution from the correction due to the larger acceptance used in 
this analysis. This uncertainty stems from the range of IP^ used in the analysis (here |IPj;| < 15cm). At 
larger absolute values of IP^ the acceptance correction for the SPD tracklets grows, and the uncertainty 
with it, being therefore T]-dependent and largest at |t]| ss 2. 

The various sources of systematic uncertainties for the FMD measurements are detailed elsewhere [6], 
but will be expanded upon in the following since some values have changed due to better understanding 
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of the detector response. 

In the analysis, three T]-dependent thresholds are used. The values for these thresholds are obtained by 
fitting a convoluted Landau-Gauss distribution [18] to the energy loss spectrum measured by the FMD 
in a given T] range. The uncertainties associated with these thresholds are detailed below. 

A charged particle traversing the FMD can deposit energy in more than one element i.e., strip, of the 
detector. Therefore it is necessary to re-combine two signals to get the single charged-particle energy 
loss in those cases. This recombination depends on a lower threshold for accepting a signal, and an upper 
threshold to consider a signal as isolated i.e., all energy is deposited in a single strip. The systematic 
uncertainties from the recombination of signals are found by varying the lower and upper threshold 
values within bounds of the energy loss fits and by simulation studies. 

To calculate the inclusive number of charged particles, a statistical approach is used [6] . The strips of the 
FMD are divided into regions, and the number of empty strips is compared to the total number of strips 
in a given region. Strips with a signal below a given threshold are considered empty. The threshold was 
varied within bounds of the energy loss fits and investigated in simulation studies to obtain the systematic 
uncertainty. 

The data-driven correction for secondary particles defined in Eq. 2 is derived from fhe previously pub¬ 
lished resulfs, and as such confains confribufions from fhe sysfemafic uncerfainfies of fhose resulfs [6]. 
Facforing ouf common correlated uncerfainfies e.g., fhe confribufion from fhe cenfralify deferminafion, 
we find a confribufion of 4.7% from fhe previously published resulfs. By sfudying fhe variafion of fhe nu- 
merafor of Eq. 1 under differenf experimenfal conditions e.g., differenf dafa-faking periods, and adding 
fhe variance in quadrafure, fhe uncorrelafed, fofal uncerfainfy on E(ri) is found fo be 6.1%. Systematic 
uncerfainfies can in general not be cancelled befween fhe numerafor and denominator of Eq. 1 , since fhe 
same T] regions are probed by differenf defecfor elemenfs in each. 

Nofe, fhaf fhe previously published resulf [6] used in Eq. 1 already carries a 2% sysfemafic uncerfainfy 
from fhe particle composition and pj disfribufion [6]. This confribufion is confained in fhe 4.7% quoted 
above, and is propagated fo fhe final 6.1% systematic uncerfainfy on E(ri). 

Finally, if was found fhrough simulafions fhaf fhe accepfance region of FMD 1 is parficularly affected 
by fhe variations in fhe number of secondary parficles sfemming from variafions in fhe particle compos- 
ifion and pj disfribufion, and gives rise fo an additional 2% sysfemafic uncerfainfy, which is added in 
quadrafure fo fhe resf of fhe sysfemafic uncerfainfies, buf only for rj > 3.7. 


Detector 

Source 

Uncertainty 

(%) 

Common 

Centrality 

0.4-6.2 

SPD 

Background subtraction 

0.1 


Particle composition 

1 


Weak decays 

1 


Extrapolation to pT = 0 

2 


Event generator 

2 


Acceptance 

0-2'^ 

FMD 

Recombination 

1 


Threshold 

+ 1 
-2 


Secondary particles 

6.1 


Particle composition & pi 

2* 


Table 2: Summary of systematic uncertainties: the common systematic uncertainties shared by both the SPD 
and the FMD, and the uncertainties particular to the detectors. ^Pseudorapidity dependent uncertainty, largest at 
|tj I =2. ^Additional contribution in 3.7 < 7] <5. See also text. 
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5 Results 

Figure 2 shows the charged-particle pseudorapidity density for different centralities from each detector 
separately. 



^7 


Pb-Pb^Aii;^ = 2.76TeV 
0-5% 

5-10% 



10 - 20 % 

20-30% 

30-40% 

40-50% 

50-60% 

60-70% 

70-80% 

80-90% 

FORWARD 


II 11111111 Uncorr.syst.unc. 

Corr.syst.unc. 
—TRACKFETS 



Uncorr.syst.unc. 

Corr.syst.unc. 

PFB 726,610 (2013) 
PRF 106,032301 (2010) 
Corr.syst.unc. 


Figure 2: [Colour online] Measurement of dAch/dp per centrality from SPD (squares) and FMD (circles) separ¬ 
ately. Error bars reflect the total uncorrelated systematic uncertainty and statistical error on each point. Error bars 
on the left and right reflect the correlated systematic uncertainties on the SPD and EMD points, respectively. Previ¬ 
ously published results for 0-30% over the full pseudorapidity range (diamonds) [6] and for 0-80% at mid-rapidity 
(stars) [7] are also shown. 

The combined distributions in Fig. 3 are calculated as the average of the individual measurements from 
the FMD and SPD, weighted by statistical errors and systematic uncertainties, omitting those which 
are common such as that from the centrality determination. The distributions are then symmetrised 
around p = 0 by taking the weighted average of ±p points. Points at 3.5 < p < 5 are reflected on to 
—5 < p < —3.5 to provide the dAch/dp distributions in a range comparable to the previously published 
results [6]. 

The lines in Fig. 3 are fits of 

— lul — t 2^ 

/GG(fl;Ai,ai,A2,a2) =Aie -A2e "''2 , (5) 

to the measured distributions. The function /gg is the difference of two Gaussian distributions centred 
at p = 0 with amplitudes Ai, A 2 , and widths ai, (72- The function describes the data well within the 
measured region with a reduced smaller than 1. We find values of A 2 /A 1 for all centralities, from 
0.20 to 0.31 but are consistent within fit uncertainties, with a constant value of 0.23 ±0.02. Fikewise 
values of a 2 /cJi for all centralities, ranges from 0.28 to 0.36 and are consistent with a constant value of 
0.31 ±0.02. 

Qualitatively the shape of the charged-particle pseudorapidity density distributions broadens only slightly 
toward more peripheral events, consistent with the above observation. Indeed, the full-width half¬ 
maximum (FWHM) shown in Fig. 4 versus the number of participating nucleons (Apart) — calculated 
using a Glauber model [13] — increase sharply only in the very most peripheral collisions. The dAch/dp 
distributions does not extend far enough to calculate reliable values for FWHM directly from the data. 
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Figure 3: [Colour online] Measurement of dA(;h/dp for all centralities and a broad rj range. Combined and 
symmetrised dA(;h/d77 over 30-90% centtality from both SPD and FMD (circles). Open boxes reflect the total 
uncorrelated systematic uncertainties and statistical errors, while the filled boxes on the right reflect the con'elated 
systematic uncertainty. Also shown, is the reflection of the 3.5 < Tj < 5 values around rj = 0 (open circles). 
Previously published results for 0-30% over the full pseudorapidity range (diamonds) [6] are also shown. The 
lines correspond to fits of Eq. 5 to the data. 


Instead fooiv) ~ ^^^(/gg )/2 = 0 was numerically solved, and the uncertainties evaluated as the error 
of /gg at the roots, divided by the slope at those roots. The width of the dAch/dt] distributions follows 
the same trend, in the region of 0-50%, as was seen in lower energy results from PHOBOS reproduced 
in Figure 4 for comparison [2]. 

Figure 5 presents the charged-particle pseudorapidity density per average number of participating nuc¬ 
leon pairs (Apart) /2 as a function of the average number of participants (Apart)- Although there is a slight 
increase in the ratio to the central pseudorapidity density distribution at low (Apart) (see lower part of 
Fig. 5), the uncertainties are large and no strong evolution of the shape of the pseudorapidity density 
distribution over pseudorapidity with respect to centrality is observed. The ratio at 3.5 < |T7| < 4.5 does 
deviate somewhat in peripheral collisions, which is attributed to the general broadening of the pseu¬ 
dorapidity density distributions in those collisions. 

To extract the total number of charged particles produced in Pb-Pb collisions at various centralities, a 
number of functions, including Eq. 5, is fitted to the dAch/dt] distributions. A trapezoid 


fy(V A X < 


® 1^1^ Tbeam 

(Tbeam % ^) V ^ ^ 

(Tbeam A/) 11] | < M 

(Tbeam-t?) r] > +M 


( 6 ) 


was successfully used by PHOBOS to describe limiting fragmentation [2]. Here, [—M,M] is the range 
in which the function is constant, and A is the amplitude. The parameterisation 


/p(77;A,a,j3,a) =A 


1 — 1/ [acosh(T])]^ 
1 _pg(|t?|-/3)/a 


(V) 
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Figure 4: [Colour online] Full-width half-maximum of the charged-particle pseudorapidity distributions versus 
the average number of participants. The uncertainties on the ALICE measurements are from the fit of /gg only 
and evaluated at 95% confidence level. Also shown are lower energy results from PHOBOS [2]. 
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Figure 5: [Colour online] The charged-particle pseudorapidity density distributions scaled by the average number 
of participants in various pseudorapidity intervals as a function of the number of participants. The four right-most 
points (open symbols) in each 77 range, as well as the mid-rapidity points (circles) are from previously published 
results [ 6 , 17]. The uncertainties on (Apart) from the Glauber calculations are only included on the points at mid¬ 
rapidity. Thus, the uncertainty band around the mid-rapidity points reflect both the measurement uncertainties and 
the uncertainty on (Apart), while other 7] ranges only show the measurement uncertainties. The lower part shows 
the ratio of each disttibution to the previously published distributions for j?] | < 0.5. 
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as suggested by PHOBOS, is likewise fitted to the dAch/dt] distributions. The parameter a expresses the 
width of the distribution, and a and jS, and expresses the width and depth of the dip at T] « 0, respectively. 
A is an overall scale parameter. Finally, to remedy some of the obvious defects of the trapezoid i.e., a 
non-continuous first derivative at T] = M, we use a Bjorken-inspired function [6] 


/B(T];A,At,a) =Ax < 


2a^ 

2ct2 


T] < -At 
T] > +At 

Ifil 


( 8 ) 


which has plateau at A for |t]| < jx connected to Gaussian fall-off beyond ±/t. The fitted functions 
are integrated over T] up to the beam rapidity iybeam = ±7.99. Although the dAch/dt] distributions in 
principle continue to infinity, there is no significant loss in generality or precision by cutting the integral 
at T] = ±ybeam sincc the distributions rapidly approach zero. Notice that all parameters of the functions 
are left free in the fitting procedure. All functions give reasonable fits (with a reduced smaller than 1), 
though the trapezoid and Bjorken-inspired ansatz are too flat at the mid-rapidity. The calculation of the 
central values and uncertainties are done as for previous results [6]: The central value is calculated from 
the integral of the trapezoid fit to compare directly to previous results; the spread between the integrals 
and the central value is evaluated to obtain the uncertainty on the total Ach- 


The extrapolated total Ach versus (Apart) is shown in Fig. 6, and compared to lower energy results from 
PHOBOS [19]. At FHC energies the particle production as a function of (Apart) shows a similar beha¬ 
viour to the lower energy results, and the factorisation [2] in centrality and energy seems to hold (see fit 
in Fig. 6). 



Figure 6: [Colour online] Exttapolation to the total number of charged particles as a function of the number of 
participating nucleons [13]. The uncertainty on the extrapolation is smaller than the size of the markers. The 
four right-most points are the previously published results [6]. A function inspired by factorisation [2] is fitted 
to the data, and the best fit yield a = 35.8 ±4.2, b = 0.22±0.05 with a reduced of 0.18. Also shown is the 
PHOBOS result at lower energy result [19] scaled to the ALICE total number of charged particles per participant 
at (Apart) = 180. 

In Fig. 7 we show comparisons of various model calculations to the measured charged-particle pseu¬ 
dorapidity density as a function of centrality. The centrality class for a given model-generated event was 
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determined by sharp cuts in the impact parameter b and a Glauber calculation [13]. 

The HIJING model [3] (version 1.383, with jet-quenching disabled, shadowing enabled, and a hard px 
cut-off of 2.3 GeV) is seen to overshoot the data for all centralities. In addition, the distributions at all 
centralities decrease with increasing | rj | faster than the data would suggest. 

AMPT [4] without string melting reproduces the data fairly well at central pseudorapidity for the most 
central events — exactly in the region it was tuned to, but it fails to describe the charged-particle pseu¬ 
dorapidity density for more peripheral events. Also, AMPT without string melting would suggest a wider 
central region than supported by data, and similarly to HIJING decreases faster than the data. AMPT 
with string melting — which essentially implements quark coalescence, and therefore a more predom¬ 
inant parton phase — is seen to be very flat at mid-rapidity and underestimates the yield, except for 
peripheral collisions. 

Finally, EPOS-LHC [5] reproduces the shape fairly well, but underestimates the data by 10 to 30%. 

6 Conclusions 

The charged-particle pseudorapidity density has been measured in Pb-Pb collisions at .^^nn = 2.76TeV 
over a broad pseudorapidity range, extending previous published results by AEICE to more peripheral 
collisions. In the mid-rapidity region the well-established tracklet procedure was used. In the forward 
regions, a new data-driven procedure to correct for the large background due to secondary particles 
was used. The results presented here are consistent with the behaviour previously seen in more cent¬ 
ral collisions and in a limited pseudorapidity range. No strong evolution of the overall shape of the 
charged-particle pseudorapidity density distributions as a function of collision centrality is observed. 
When normalised to the number of participating nucleons in the collision, the centrality evolution is 
small over the pseudorapidity range. Since the measurement was performed over a large pseudorapidity 
range (—3.5 < rj < 5), it allows for an estimate of the total number of charged particles produced in 
Pb-Pb collisions at = 2.76TeV. The total charged-particle multiplicity is found to scale approx¬ 
imately with the number of participating nucleons. This would suggest that hard contributions to the total 
charged-particle multiplicity are small. From peripheral to central collisions we observe an increase of 
two orders of magnitude in the number of produced charge particles. A comparison of the data to the 
different available predictions from HIJING, AMPT, and EPOS-EHC show that none of these models 
captures both the shape and level of the measured distributions. AMPT however comes close in limited 
ranges of centrality. The exact centrality ranges that AMPT describes depend strongly on whether string 
melting is used in the model or not. EPOS-EHC — although systematically low — shows a reasonable 
agreement with the shape of the measured charged-particle pseudorapidity density distribution over a 
wider pseudorapidity range. 
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